arXiv: 1505.04747v2 [cond-mat.quant-gas] 17 Jun 2015 


Cooling of a one-dimensional Bose gas 
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We experimentally study the dynamics of a degenerate one-dimensional Bose gas that is subject 
to a continuous outcoupling of atoms. Although standard evaporative cooling is rendered ineffective 
by the absence of thermalizing collisions in this system, we observe substantial cooling. This cooling 
proceeds through homogeneous particle dissipation and many-body dephasing, enabling the prepa¬ 
ration of otherwise unexpectedly low temperatures. Our observations establish a scaling relation 
between temperature and particle number, and provide insights into equilibration in the quantum 
world. 
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Introduction .—Long coherence times, the tunability of 
many parameters and the ability to precisely probe and 
manipulate their quantum state make ultracold atomic 
gases a very promising and versatile tool to study the 
physics of quantum many-body systems [1]. The stan¬ 
dard technique to reach the necessary ultracold temper¬ 
atures in these systems is evaporative cooling nm, a fun¬ 
damental process present in many physical systems from 
hot liquids to stellar clusters. It relies on the selective 
removal of the most energetic particles from a trapped 
gas and the subsequent re-thermalization of the gas to a 
lower temperature through elastic collisions. For efficient 
cooling this cycle is repeated continuously, increasing the 
phase-space density of a gas at the price of reducing the 
total number of atoms. 

Intuitively, systems that do not thermalize can not 
be cooled through particle dissipation. In this letter we 
demonstrate that this notion is incomplete by studying 
the dissipative dynamics of a one-dimensional (ID) Bose 
gas with contact interactions. We observe a substantial 
decrease in temperature as a result of homogeneous par¬ 
ticle dissipation, reaching temperatures far below huj±, 
the energy characterizing the transverse confinement. 
In this deep ID regime thermalization is strongly sup¬ 
pressed W- Thermalizing two-body collision are frozen 
out [7] and three-body collisions mi or phonon-phonon 
scattering can be neglected on our timescales. We 

find that the observed cooling can be modeled as a con¬ 
tinuous density reduction extracting energy from the den¬ 
sity quadrature of the free phononic excitations. To¬ 
gether with a continuous many-body dephasing this re¬ 
duces the occupation number of each phonon mode and 
leads to a colder system. 

Experiment .—Our experimental system is a ID Bose 
gas of ®^Rb atoms prepared in the anisotropic magnetic 
trapping potential of an atom chip [12) . Initially, a pre¬ 
cooled three-dimensional (3D) cloud of thermal atoms, 
prepared in the |F, mi?) = \2,2) ground state is loaded 


into the trapping potential. The gas is then cooled 
through the condensate transition and into the ID regime 
by conventional evaporative cooling. The evaporation 
is realized through a weak driving of energy-selective 
radio-frequency (RF) transitions to untrapped Zeeman 
states m- The cooling ramp of the RF frequency is 
performed over 1.6 s, with the final part being particu¬ 
larly slow to minimize collective excitations in the gas. 
The ID regime is reached when both the temperature 
and the chemical potential drop below huj±, the energy 
spacing between the ground state and the first excited 
state of the transverse confinement. At this point the 
condensate typically consists of 7000 to 10000 atoms at 
an initial temperature T of 30 to 100 nK. The result¬ 
ing macroscopic quantum system is a quasi-condensate 
described by a macroscopic wave function with a fluctu¬ 
ating phase [13]. 

This trapped ID quasi-condensate is the starting point 
of our experiments. To investigate the effects of particle 
dissipation we monitor the further evolution of the sys¬ 
tem under continuous driving of transitions to untrapped 
states. The outcoupling rate is tuned in such a way that 
atoms are removed slowly from the quantum degener¬ 
ate gas to avoid collective excitations. In contrast to 
conventional evaporative cooling the atoms are not out- 
coupled energy-selectively but at a nearly homogeneous 
rate. This is experimentally demonstrated in Fig. ^ by 
extracting a visible fraction of the gas using a single out¬ 
coupling pulse and analyzing its relation to the source 
cloud. The outcoupled atoms show the same average den¬ 
sity profile as well as the same density speckle patterns, 
indicating a coherent and homogeneous outcoupling pro¬ 
cess. 

Information about the system is extracted through ab¬ 
sorption imaging in time-of-flight, transversal to the ID 
axis |15j . To improve the individual images we employ 
a fringe removal algorithm |16j . The imaging provides 
access to the evolution of the atom number, the den- 



2 



FIG. 1. (color online). Analysis of the outcoupling mech¬ 
anism. (a) A typical density speckle pattern formed by a 
ID quasi-condensate expanding for 10.5 ms in time-of-flight. 
2.5 ms before the gas was released a small fraction was out- 
coupled by a single 2 ms RF pulse. The outcoupled cloud is 
visible below the source cloud, showing the same speckle pat¬ 
tern. (b) Comparing the longitudinal density profiles of the 
source cloud (blue) and the outcoupled cloud (red) averaged 
over many realizations reveals the same form of the profile, 
indicating a homogeneous outcoupling process, (c) Sketch 
of the outcoupling mechanism. Transversal to the ID axis 
the atoms occupy the ground state of the harmonic potential 
(blue). Outcoupling is achived by coupling these atoms to 
untrapped states, e.g. |2,0) (red) which only feels the lin¬ 
ear gravitational potential. The rate of this transition, taking 
only the transversal degrees of freedom into account, is plotted 
in (d) over the detuning of the driving field, given in units of 
the longitudinal trap frequency. The large width of this rate 
compared to the energy scale of the longitudinal dynamics 
leads to a nearly homogeneous outcoupling of atoms. 


sity profile and the temperature of the gas. The lat¬ 
ter is extracted from the longitudinal density speckle 
patterns forming in the time-of-flight expansion seen in 
Fig.[^ HZIIIH]. These patterns are a direct result of the 
longitudinal phase fluctuations in the trapped gas. Ex¬ 
tracting their normalized auto-correlation function and 
comparing it to simulated data generated through a 
stochastic Ornstein-Uhlenbeck process [19] allows us to 
infer the temperature of the gas. 

Our results for the evolution of the atom number and 
the temperature for different initial conditions, outcou- 



FIG. 2. (color online). Three different examples for mea¬ 
surements of temperature T over atom number N-. The green 
triangles are obtained through a continuous RF-induced out¬ 
coupling from a harmonic trap with trap frequencies uj± = 
2 tvx (2.1 ±0.1) kHz and wn =27rx(ll±l) Hz over a timescale 
of 25 ms. The blue circles differ in that they start from a much 
colder initial temperature and were measured over a timescale 
of 90 ms in a transversally slightly anharmonic trap created 
through a dressed state potential |20| . There the trap frequen¬ 
cies are uj± = 2n x (1.6 ± 0.1) kHz and wn = 27r x (8 ± 1) Hz. 
Finally, the red diamonds were acquired through a very differ¬ 
ent measurement procedure, where MW radiation was used 
to couple atoms from the dressed trap directly to the anti¬ 
trapped 11,1) state. Furthermore, in contrast to the other 
data points the temperature was measured after a constant 
period of outcoupling with a different outcoupling rate for 
each data point. The temperature error bars give the 68% 
confidence interval obtained from a bootstrap. The inset de¬ 
picts the speckle correlations of the coldest red and blue data 
points together with their respective thermal fits still showing 
good agreement. 


pling rates and measurement procedures are shown in 
Fig.i The outcoupling is achieved through RF tran¬ 
sitions within the F = 2 manifold or microwave (MW) 
transitions to the anti-trapped |1,1) state. We find that 
in all cases the correlation functions of the density speck¬ 
les in time-of-flight remain close to their thermal form for 
all data points. This is remarkable, as the system is close 
to an integrable point and thus can not thermalize on ex¬ 
perimentally relevant timescales. In addition, the nearly 
homogeneous, energy independent out-coupling process 
would intuitively not lead to cooling. Nevertheless, we 
observe a significant decrease in temperature down to 
ksT ^ 0.1fiuj±. Moreover, even though the system con¬ 
stantly loses atoms the ratio kgT/fj, between the thermal 
energy and the interaction energy given by the chemical 
potential drops as well. In the coldest measurement 
presented in Fig. (blue circles) it reaches values as low 
as ~ 0.25. The thermal coherence length At is inversely 
proportional to this ratio showing that the system be- 
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FIG. 3. (color online). The data from Fig.j^is plotted rescaled 
to their respective initial temperature To and atom number 
No values revealing a clear linear scaling. 

comes more coherent under dissipation. 

Fig. [^suggests that the temperature decreases linearly 
with the number of atoms. Rescaling the data points 
to their respective initial values Tq and iVo collapses the 
measurements to a single line, as shown in Fig. This 
suggests a linear scaling relation between temperature 
and atom number 


which is particularly interesting as these measurements 
are obtained from different experimental procedures and 
initial conditions. 

Model .—In a pure ID setting two-body collisions do 
not lead to a re-distribution of energy and momentum. 
However, in a quasi-lD trap with a tight transverse con¬ 
finement this condition can be broken in collisions where 
there is enough energy available to access transverse ex¬ 
cited states. These thermalizing two-body collisions are 
suppressed by a factor exp(—2 ?iu;j_/A:bT) [7] in a gas of 
non-degenerate bosons. For a degenerate gas the Bose en¬ 
hancement of the ground state leads to an even stronger 
suppression. Consequently, these collisions freeze out as 
soon as the gas enters the ID regime. Other processes 
that can lead to thermalization in our system are three- 
body collisions 1311] or phonon-phonon scattering [SHn]. 
However, their expected thermalization time-scales are 
beyond the times probed in our experiment and can not 
explain the observed cooling. A different case of evapo¬ 
rative cooling in a ID harmonic trap has been modeled 
by Witkowska et al. EH, assuming that atoms are out- 
coupled energy-selectively at the two ends of the cigar¬ 
shaped trap, which is not the physical situation in our 
experiment (see Fig. [^. 

To find a simple mechanism consistent with our exper¬ 
imental observations we start by considering the outcou- 
pling process. The rate of outcoupling via RF or MW 
transitions is determined by the overlap of the conden- 




FIG. 4. (color online). Dissipative cooling mechanism visual¬ 
ized through the Wigner function (a) of the phase and den¬ 
sity quadrature (j)j and Hj of a single mode j, which rotates 
with the eigen frequency ojj. The homogeneous outcoupling 
of atoms squeezes the Wigner function in the density quadra¬ 
ture, indicated by the black arrows, thereby reducing the en¬ 
ergy of the mode. The rotation with Dj together with a slow 
continuous squeezing leads to a (nearly) symmetric narrowing 
of the Wigner function (b) reflecting a reduced occupation of 
the mode. In our multimode system this happens in all modes 
at the same time, leading to a reduced temperature. 


sate wavefunction with the free particle states of the un¬ 
trapped Zeeman levels (see Fig. &)• The observed ho¬ 
mogeneity can be explained by the fact that the energy 
of the longitudinal degrees of freedom is not resolved by 
this process since the timescale on which atoms leave the 
trap is determined by the much tighter transversal con¬ 
finement (see Fig. W)- This means that each atom has 
the same probability to leave the trap irrespective of its 
position or energy. 

The dynamics resulting from such a homogeneous par¬ 
ticle dissipation can be intuitively understood within a 
Luttinger-liquid picture [22j . describing the low energy 
dynamics of the underlying Lieb-Lininger model [23] . 
The elementary excitations in this model are non¬ 
interacting phonon modes. Each of these modes con¬ 
tributes to the fluctuations in the gas through a density 
and a phase quadrature, in close analogy to the posi¬ 
tion and momentum quadratures of a harmonic oscilla¬ 
tor. The free evolution of such a mode j with energy 
hujj can be visualized as a rotation of the corresponding 
Wigner function with the frequency ojj. In this picture, a 
sudden homogeneous outcoupling of atoms leads not only 
to a decrease in average density but also to a correspond¬ 
ing decrease in the density fluctuations. Such an instan¬ 
taneous density reduction therefore extracts energy from 
the density quadrature of the phonon modes while leav¬ 
ing the phase quadrature unchanged. The system reacts 
to this reduction by dephasing, redistributing the remain¬ 
ing energy between the quadratures. A related process 
has recently been studied in detail using matter-wave in¬ 
terferometry after coherently splitting a ID Bose gas and 
has been shown to lead to prethermalization PI21] and 
light-cone-like spreading of thermal correlations [IS] . 
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The experimentally investigated situation of a slow, 
continuous dissipation can be modeled as a series of small 
density reductions each followed by an immediate de¬ 
phasing. This leads to a continuous decrease of the quasi¬ 
particle mode occupation, as illustrated in Fig. The 
mode occupations reacts to a sudden decrease in density 
p —t p' as n'j = + 1) rij, with a = p'/p being the 

ratio by which the average density is reduced. For a series 
of small reductions under continuous dephasing we can 
expand this factor around a = 1, which, to leading order, 
gives n'- = rij . In the case of a harmonically trapped 

gas in the Thomas-Fermi limit where the central density 
scales as p'/p = this leads to a linear scaling 

of rij with the particle number N. As the degenerate 
system is still sufficiently hot, the observable low-energy 
modes are dominated by the Rayleigh-Jeans part of the 
Bose-Einstein distribution where the mode occupation 
is proportional to the temperature rij ~ ksT/cj, with 
tj being the mode energy. Therefore, a homogeneous 
dissipation ultimately leads to a linear scaling relation 
between temperature and atom number T/Tq = N/Nq, 
as observed in Fig. An extensive theoretical analysis 
of this process can be found in |2fij . Inhomogeneities in 
the outcoupling rate, stemming for example from the en¬ 
ergy of the longitudinal degrees of freedom would lead 
to stable non-thermal energy distributions described by 
generalized Gibbs ensembles miiiH]- 

This result extends scaling laws for the dynamical com¬ 
pression and expansion of a trapped ID Bose gas |53] to 
the realm of open quantum systems. Its universality is 
a result of the phononic nature of the excitations. Sup¬ 
porting our simple model, the scaling law of Eq. 0 can 
also be found by numerically solving the Gross-Pitaevskii 
equation with an additional dissipative term. Furthere- 
more, it is interesting to note that this cooling mechanism 
is similar to cooling by adiabatic expansion, in that the 
entropy per particle stays constant. 

It has been suggested that vacuum fluctuations enter 
the system through the dissipation process US 1301 EH, 
which would lead to a reduced cooling efficiency. How¬ 
ever, all our measurements are consistent with the scaling 
of Eq. 0 derived from our simple model. 

Conclusion .—We have experimentally observed the 
cooling of a degenerate atomic gas in the ID regime by 
contious outcoupling of particles. In this process we reach 
temperatures as low as ~ 0.25 /r and ~ 0.1 huj± 
which is far below the region where thermalizing two 
body collisions freeze out, rendering standard evapora¬ 
tive cooling ineffective. The observed cooling allows us 
to go deep into the ID regime and is therefore of direct 
practical importance for all experiments with ID Bose 
gases. Our model also explains how previous experiments 
could observe such low temperatures |32l |33] . 

Finally, our simple model suggests that the mecha¬ 
nism behind the observed cooling is not limited to ID 
and could be relevant in 2D and 3D settings at very low 


temperatures where atoms are homogeneously outcou- 
pled from the quantum degenerate gas and the dephasing 
of excitations is faster then their thermalization. 
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